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The exact nature of the chiral crossover in QCD is still under investigation. In V/ = 2 and 
Nf = (2+1) lattice simulations with staggered fermions the expected 0(A/")-scaling behavior was 
observed. However, it is still not clear whether this behavior falls into the 0(2) or 0(4) universality 
class. To resolve this issue, a careful scaling and finite-size scaling analysis of the lattice results are 
needed. We use a functional renormalization group to perform a new investigation of the finite-size 
scaling regions in 0(2)- and 0(4)-models. We also investigate the behavior of the critical fluctuations 
by means of the 4^^-order Binder cumulant. The finite-size analysis of this quantity provides an 
additional way for determining the universality class of the chiral transition in lattice QCD. 


I. INTRODUCTION 

Quantum Chromodynamics (QCD) at finite tempera¬ 
tures and chemical potential is currently one of the most 
actively researched topics in theoretical physics. In par¬ 
ticular, the exact nature of the transition from hadronic 
phase to quark-gluon plasma is of great importance for 
interpretation of experimental results from heavy ion col¬ 
lisions [1] and for the understanding of the evolution of 
the universe in its early stage. 

Two transitions take place in QCD: the confinement- 
deconfinement and the chiral transition. The first one 
turns out to be a first order phase transition for pure 
gluonic systems, while it becomes a crossover in the pres¬ 
ence of quarks. In contrast, the nature of the chiral tran¬ 
sition depends on the number of dynamical quark fla¬ 
vors, but also on the strength of the explicit symmetry 
breaking and the strength of the chiral anomaly [HIS]. 
For zero chemical potential, vanishing anomaly and two 
massless quark flavors it is expected to be a phase tran¬ 
sition of second order. In this case QCD falls into the 
0(4) universality class [2]. For massive quarks, however, 
the phase transition disappears and we expect to observe 
a crossover which should still fulfill the 0(4) scaling be¬ 
havior. In the presence of the strong chiral anomaly and 
for two quark flavors in the chiral limit one expects that 
the chiral phase transition is of first order ElISl. 

A convenient method to study full QCD are lattice sim¬ 
ulations. One particular way to implement fermions on 
the lattice is to use staggered fermions. This method al¬ 
lows to reproduce pion masses which are close to the chi¬ 
ral limit. However, staggered fermions exhibit only a re¬ 
duced chiral symmetry and we can also expect to observe 
a chiral transition governed by 0(2)-symmetry. Finite 
simulation volumes poses an additional problem: contin¬ 
uous chiral symmetry cannot be broken spontaneously in 
a finite volume and explicit symmetry breaking in form 
of finite quark masses is mandatory. So, no real continu¬ 
ous phase transition can be observed in such simulations 
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but only a crossover. These facts lead to complications in 
the interpretation of lattice data and hence, to ambigu¬ 
ous results: For a long time, data provided by lattice 
simulations with two dynamical or (2+1) quarks using 
the staggered formulation did not exhibit the expected 
©(A/")-scaling [4]-[6] or gave evidence for a first-order phase 
transition mm- Results from present-day calculations 
are in a very good agreement with ©(A/")-scaling behavior 
[MI]. Nonetheless, there is still an ambiguity about the 
exact universality class of the transition observed. How¬ 
ever, especially the finite simulation volumes can affect 
the universal behavior of the thermodynamic observables 
in a very specific way, depending on the nature of the 
transition in the limit V ^ oo. Therefore, finite-volume 
effects could help to shed light on the nature of the chiral 
crossover observed in the lattice QCD [91[I2|. 

The chiral transition seems to be a continuous one. In 
this case long-range fluctuations are dominant and the 
microscopic details of the system are no longer relevant. 
The observed behavior is then basically determined by 
symmetries and dimensionality or, with other words, by 
the universality class of the system. It means that close to 
the transition point the thermodynamic quantities obey 
power-laws and scaling functions which are characterized 
by universal critical exponents. Even in the presence of 
finite volumes the behavior of the thermodynamic observ¬ 
ables can be described by still universal finite-size scaling 
functions. 

Since long-range fluctuations play such a prominent 
role at continuous transitions, we can try to investigate 
the nature of the chiral transition in lattice QCD us¬ 
ing the higher-order fluctuations. For this purpose we 
can use a quantity called Binder cumulant [T3|. The 
Binder cumulant has already been used successfully for 
determining the critical value of the quark mass for three 
degenerated massive flavors, where the phase transition 
is continuous and falls into the Z(2) universality class 

HMg. 

So far, scaling functions, finite-size scaling functions 
and Binder cumulant have been determined mainly by us¬ 
ing ©(A/")-symmetric spin-model lattice simulations [16]- 
|20]. These results were often used in the analysis of 
lattice QCD data [6] [H [161 [2TII23] . In the present pa- 
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per we use an alternative technical framework, the func¬ 
tional renormalization group approach (FRG). We cal¬ 
culate scaling and finite-size scaling functions for the 3- 
dimensional linear 0(2) model and compare our results 
with findings from [24] and m where the 0(4) model was 
considered. For both models we also calculate the Binder 
cumulant and investigate its behavior in finite volumes. 
Our approach is complementing 0(A/')-spin-model lattice 
simulations. Our results can be directly applied in the 
scaling analysis of the lattice QCD data. 

This paper is organized as follows: In Sec. we briefly 
discuss the setup of our FRG formalism and introduce 


scaling and finite-size scaling functions. Then in Sec. Ill 


we define Binder cumulant of the 4^^ order and discuss 
its properties. After, in Sec. W and Sec. [V| we present 
our results for critical exponents and scaling functions 
in infinite and finite volumes. In Sec. EH we show and 
discuss numerical results for the Binder cumulant. Our 
concluding remarks can be found in Sec. |vnl 


II. METHOD AND MODEL 


A. FRG applied to 0(A)-models 


In our investigations we use the functional renormaliza¬ 
tion group applied to continuous 0(A^)-symmetric mod¬ 
els in 3 spatial dimensions. At the ultra-violate (UV) 
scale A the bare action of these models is defined as 

T4cf>\ = Id^x{^Z4d^cf>f + UA{cf> ')) , ( 1 ) 


where 0 = (a, tti, ..., represents multiple, 

effective scalar degrees of freedom, is a wave function 
renormalization and an effective potential de¬ 

fined at the scale A. This potential should depend only 
on the powers of and should exhibit 0(A/')-symmetry. 
If we introduce the finite explicit symmetry breaking in 
our calculations, we add the term —Her to the poten¬ 
tial. In our calculations we use the so-called local poten¬ 
tial approximation (LPA) where we assume that expec¬ 
tation values of the fields do not have any spatial depen¬ 
dence and the wave function renormalization is constant: 

= 1. This assumption leads to vanishing anomalous 
dimension, = 0. Even though the value of r] for 0{N)- 
models was measured to be finite, it is still relatively 
small compared to one, see, e.g.. Ref. [25]. Therefore our 
assumption of constant wave function renormalization is 
well justified. 

The RG flow of the effective action can be described 
by the Wetterich flow equation [26] : 


dtrk = 




( 2 ) 


fulfill certain constraints [5S]. We are free in the choice 
of R/e, so we can use it for the optimization of the flow 
PTHSO] . In this work we use the 3-dimensional Litim’s 
optimized regulator [29] given by 

Rk{p^) = {k^-p'^)0{k^-p'^) . (3) 


Using this regulator function and applying the LPA we 
find the following flow equation oiUk for the case of in¬ 
finitely large volumes: 


dtUk[cp^] 


r 




( 4 ) 


The scale-dependent expressions Mj and are defined 
as the eigenvalues of the second-derivative matrix of Uk 
and are equal to the square of bosonic masses in the limit 
k ^ 0. Note that they still depend on the field 0. 

Since we do not know the exact form of the effective 
potential U/c, we expand it around its particular mini¬ 
mum (j) = (^o(k), 0,..., 0)^: 


M 

Uk=J2a^(k)(cP^-a^(k)r-H-a, (5) 

m=0 

where cro(k) ^ 0 only if the symmetry is broken. If we 
also introduce a finite explicit symmetry-breaking, we 
choose 7^ 0, otherwise H = f) 

Expansion of the right-hand side of RG-flow, Eq. 
around the minimum of Uk will provide us with a set of 
highly coupled differential equations for couplings am{k). 
In order to fix the scale-dependent expectation value of 
field <j, we use an additional condition which ensures that 
we are expanding our potential around the actual physi¬ 
cal minimum: 


dUk 

da 


= 0 

cr=ao(k),7r^=0 


2ai{k)ao{k) = H . (6) 


Now we have all information we need to solve the RG 
flow. However, the number of allowed couplings am{k) 
is infinite. So we have to choose an appropriate trun¬ 
cation scheme for our potential. Unfortunately there is 
no argument allowing us to neglect couplings of higher 
order a priori like it can be done in perturbative calcu¬ 
lations. However, in Eais] it was found that inclusion 
of a small number of couplings is already sufficient for 
0{N)-modeh. In our present calculations we have used 
M = 6 and have checked explicitly that our results do 
not change considerably if we use additional couplings. 

At the Wetterich flow equation for finite volume, we 
replace spatial momenta integrations by a sum over dis¬ 
crete momenta. In this step it is possible to use periodic 
or anti-periodic boundary condition for discretized mo¬ 
menta. Since most of lattice QGD simulations use peri¬ 
odic boundary condition, we also use this: 


where t = \n{k/A). The function Rk in the expression 
above is the so-called regulator function which controls 
the Wilsonian momentum-shell integrations and has to 



— oo 


E 

rii 


with i = 1, 2, 3 (7) 
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with 


P 


2 


47r^ 

IX 



(8) 


The flow of our effective potential 0 is then modified to 


dtUk[<l)^] = k- 


1 


+ 


(N-l) 


k^ + M^ k^ + M^ 


B{kL) , (9) 


where B{kL) is a mode-counting function which includes 
all information we need for fluctuation modes allowed in 
a particular volume with extent L. It is provided by 

BikL) = - P • ( 10 ) 

^ ^ nez^ 


The limiting behavior of this function is important for 
understanding the contributions of different fluctuation 
modes to the RG-flow: Since we use periodic boundary 
condition in our calculations, the mode counting function 
behaves in the limit kL ^ 0 as: 


lim B(kL) 

kL^O 


1 

(kL)^ ‘ 


( 11 ) 


This is due to existence of zero-momentum mode for the 
choice of periodic boundary condition. For very small 
volumes the dynamics of the system are basically gov¬ 
erned by this zero-momentum mode. 

For infinitely large volumes the sum in Eq. (10) be¬ 
comes the volume of a sphere with the radius r = 
kLj(2i{) and we find 


lim B(kL) = 

kL^oo 


1 

Gtt^ 


( 12 ) 


Thus, we recover the flow equation for infinite volume. In 
our finite-volume calculations we use the same expansion 
of the effective potential as given in Eq. (§. Since for 
finite volumes a finite explicit symmetry-breaking term 
is mandatory, we have to choose H 0. 

The ERG approach described above allows us to calcu¬ 
late scaling and finite-size scaling functions as well as the 
Binder cumulant. These functions are introduced below. 


B. Scaling functions in infinite volume 

The behavior in the vicinity of the critical transition 
point is governed by the free-energy density. This quan¬ 
tity consists of a singular and a regular part 

f = MT,H)+fr{T,H), (13) 

where only the singular part is responsible for critical 
behavior, whereas the regular one leads at most to some 
finite corrections. 

In our model we assume that only temperature T and 
explicit symmetry breaking H are the relevant couplings. 


In order to remove all system-specific scales, we introduce 
the rescaled temperature t and the rescaled symmetry¬ 
breaking field h as new variables: 


t = 


T-Tc 

To 


( 14 ) 

TLq 


where Tc is the critical temperature and Tq and Hq are 
system-dependent normalization constants. Then, the 
singular part of the free-energy density is a function of t 
and h: fs{T,H) = fs{t,h). 

Close to the critical point, correlation length has in¬ 
finite range ^ oo) and the system becomes scale in¬ 
variant. Therefore the singular part of the free energy 
density is invariant under a rescaling of the length with 
a factor a: 


fs{t,h)=a . (15) 


The critical exponents can be then expressed in terms of 
yt and yn' 


1 P5 

yt = - ^ yh = — • 

U V 


(16) 


Since we have two relevant couplings in our model, only 
two critical exponents are independent. All others can 
be obtained using the following scaling laws: 


7 =( 2 -r?)iy, 7 = ^(<5 - 1 ) , P - 2 + r])v , 


ud — /3(1 + (5) 


d + 2 — 77 
d — 2 T] 


(17) 


The critical exponents (3 and S govern the behavior of 
the order parameter M, which is associated in our model 
with the expectation value of the field a: 

M{t,h = 0) = {-tf , M{t = 0,h) = h^/^ . 

(18) 

The critical exponents u and 7 describe the behavior of 
the correlation length ^ and longitudinal susceptibility y 
correspondingly: 


^cx|^r^ xocllpv (19) 

In this paper we investigate only the behavior of the 
order parameter M. Using the rescaled form of the free 
energy density, Eq. M can be derived from its ther¬ 
modynamic definition: 

M = ao =-^ = h'^^^fMiz) , ( 20 ) 

where 2 : = is a new single scaling variable and 

is invariant under the rescaling of t and h. 

In the expression above, /m(^) is a scaling function for 
the order parameter M. It turns out that this function 
is universal for an universality class. Eor other thermo¬ 
dynamic observables the corresponding scaling functions 
can be derived at the similar manner. 
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C. Scaling functions in finite volume 


If a thermodynamic system is put into a finite volume, 
the correlation length is bounded from above by the sys¬ 
tem extent L. Since the critical point can be reached only 
in the limit L ^ oo, the system extent becomes an addi¬ 
tional coupling, and the critical behavior changes. Since 
the singular part of the free-energy density is now a func¬ 
tion of three variables, we need to introduce two scaling 
variables in order to define a finite-size scaling function. 

According to Fisher’s finite-size scaling hypothesis [32] . 
the ratio of thermodynamic quantities in the finite-size 
system to those in an infinite system depends on only 
the ratio ^/L. This implies that in the absence of h the 
system extent has to scale with t exactly in the same way 
as Eq. (IT^. We also know from definition of 2 ; that 
t oc These observations allow us to introduce a 

new scaling variable h* = ^ where / = L/Lq is the 

renormalized system extent, and Lq is system specific. 
Then, using Eq. ( [2Q| ) for the order parameter, we find 

M{t,h,l)=l-^/^QM{zX) , ( 21 ) 

where Q{z,h*) is a finite-size scaling function in lead¬ 
ing order. Eor completeness, we also specify the leading- 
order finite-size scaling correction |32) : 




Qm{z, h*) + ^*) + • • • 


( 22 ) 

where uj is the critical exponent associated with the first 
irrelevant operator in the renormalization group flow. 
This additional term influences the behavior of the sys¬ 
tem in the vicinity of the critical point for small L and 
needs to be removed in order to isolate the universal 
finite-size scaling function. In fact, we observe this cor¬ 
rection in our calculations for small volumes. However, 
since we use very large volumes to fix Q{z^h*)^ we will 
neglect finite-size scaling corrections and use the above 
leading-order expression for the order parameter. 


III. BINDER CUMULANT 

In addition to the universal scaling function for the 
order parameter we also investigate a higher-order fluc¬ 
tuation quantity, the so-called Binder cumulant of the 

4th 

m- In the infinite-volume limit, this quan¬ 
tity exhibits a value at the critical point which is specific 
for some particular universality class. Therefore, it is 
often used for the localization of the critical point, e.g. 

nmisiiss]. 

The universal values of the Binder cumulant for the 
three-dimensional 0(2)- and 0(4)-models have already 
been determined at high accuracy by using spin-model 
lattice simulations [191 [20]. To our knowledge, the Binder 
cumulant has never been calculated previously using 
ERG. Inj^we present such a calculation for the 0(2) and 


0(4) universality classes in LPA. Using our approach we 
can also investigate the influence of finite-volume effects 
on the Binder cumulant. 

The convenient definition of the 4^^-order Binder cu¬ 
mulant is given by 


B^ = 


((0T) 

(</>T 


(23) 


with (j) = (cr, TTi)^ and i number of Goldstone modes. 
In this quantity we compare the total contribution of all 
possible fluctuations of the 4^^ order to the contribution 
from trivial Gaussian fluctuations. 

In the phase with broken symmetry, the system ob¬ 
tains a very high stability. In this regime fluctuations 
are suppressed by powers of ^ (see and we expect 

((0^)^)^M^, (24) 


for T <^Tc. Therefore, for any 0(A)-model the value of 
H 4 should approach 1 for decreasing temperatures. 

In contrast, for T ^ Tc fluctuations become the lead¬ 
ing contributions. At large T the contribution from fluc¬ 
tuations of some particular order goes into saturation 
and H 4 approaches some finite value. Eor models with 
a different number of degrees of freedom, the number of 
fluctuations contributing to ( 0 ^) is obviously not iden¬ 
tical and the Binder cumulant H 4 approaches a value, 
which is specific for a particular model. Our calculations 
suggest a general expression for this limiting behavior for 
0 (A)-symmetric models: 


H4 = 


A + 2 
N 


(25) 


So, for T Tc we expect B 4 = 2 for 0(2)- and B 4 = 
3/2 for 0(4)-model. These values are confirmed by our 
numerical calculations. 

Additionally the Binder cumulant is directly a finite- 
size-scaling function. If we consider almost vanishing 
symmetry-breaking field this function is given by 

B4 = QB{tL^'\L-'^,...) . (26) 


Qb depends on the scaling variable and possibly on 
some other irrelevant operators which we specify only up 
to the leading order, i.e., we assume only the finite-size 
corrections proportional to L~^^ with cj > 0. Expanding 
this finite-size scaling function to the lowest obtainable 
order in both variables 


QB(tL^'\ i““) = ao + + a2L-^ + • • • • (27) 

Erom this expression we see that exactly at the critical 
point (t = 0, 1/ ^ oo) the Binder cumulant is simply 
given by a constant. This constant is in general different 
for different universality classes. Therefore it can be po¬ 
tentially applied in order to determine universality class 
of a particular system, e.g., in lattice QCD simulations. 
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0 ( 2 ) 

H = 0 


J7 = 1.0 X 10“^® 




derivative 

fit 

derivative 

fit 

final 

(3 

0.3537 

0.3536 

0.3540 

0.3538 

0.3538(2) 

7 

1.4143 

1.4139 

1.4134 

1.4142 

1.4140(4) 


0 ( 2 ) 

derivative 

fit 

7//? + 1 

(d + 2 - 77)/(d - 2 + ? 7 ) 

final 

(5 

4.9997 

4.9995 

4.9966 

5.0000 

4.9990(15) 


Table L Our results for the critical exponents /3, 7 and S for the 0(2)-model in LPA. Using different fit techniques allows us 
to estimate the uncertainty of our determinations. In the evaluation of S we also use our results for /3 and 7 and theoretical 
prediction for S. Since in our formalism 77 = 0 and d = 3, the theoretical prediction for S is 5.0. 


IV. INFINITE-VOLUME SCALING 
A. Critical exponents 

To this day, critical exponents for three-dimensional 
0 ( 2 )-model were calculated at very high accuracy us¬ 
ing different methods such as lattice Monte-Carlo sim¬ 
ulations of spins [m HD, perturbative field-theoretical 
[sms], and RG calculations EH ED [361ET]. Since we 
use LPA in this particular work which implies the van¬ 
ishing anomalous dimension 77 = 0 , we do not aim to add 
it to this list. However, in order to provide a consistent 
evaluation of scaling functions, we need to calculate crit¬ 
ical exponents within the LPA and to use them in the 
following analysis. 

In our formalism we use d = 3, therefore, we cannot de¬ 
fine temperature in the field-theoretical sense. Nonethe¬ 
less, we can find a parameter in the RG-fiow which does 
the same job as T, i.e., controls the phase transition. In 
our case it is the initial value of the expectation value of 
(j) at the UV scale A. So, we suppose the existence of an 
expansion (0o(A) - - (T - Tc) 

In our calculations we use A = 10 GeV. For this setup 
we obtain the following critical UV-value: 

^criticalA) = 37.996488987996596 MeV^/^ . (28) 

As a point of fact, this accuracy in is necessary 

in order to observe the fixed point in the RG-fiow. For 
simplicity we will use in the following the notation T and 
Tc for 00 (A) and 

We calculate the critical exponents (3 and 7 in the 
chiral limit (id = 0 ) using the power-laws discussed in 
Sec. 11111 We use two different fitting techniques: direct 
fitting and taking the numerical derivative of log(M) and 
log(x) in the limit t > 0. In order to estimate error, 
we repeat these calculations for small explicit symmetry 
breaking H = 10“^^ MeV^^^. Our results can be found 
in the upper part of Table [I[ In the following analysis we 
use averaged values of the critical exponents. 

For the determination of 5 we evaluate the order pa¬ 
rameter M exactly at the critical temperature for dif¬ 
ferent symmetry-breaking fields H. Once again, we pro¬ 


vide a direct fit of our results and a fit to the numerical 
derivative of log(M). In addition, we calculate 6 using 
the scaling law ^ = 7 /^ + 1 and our results for 0 and 7 
from Table [Tj In the determination of the averaged value 
of (5, we also use the theoretical prediction: 


d T 2 — Tj 
d — 2 + 77 


(29) 


In d = 3 dimensions and in the absence of the anomalous 
dimension (LPA) we expect 6 = b. Our critical exponents 
are summarized in the lower part of Table |Tj 

Our results seem to be consistent and agree within 
0 . 1 % with the results from m, where exact renormal¬ 
ization group in LPA was used to calculate the critical 
exponent u = 2/3 (in LPA). Never the less, we should 
keep in mind that all our calculations include some ad¬ 
ditional systematic errors corresponding to truncation of 
the effective potential U ( 0 ^) and to neglecting the kinetic 
terms of higher order in the effective action P [0] . How¬ 
ever, we have considered a relatively large number of n- 
point couplings {rimax = 12) and results in Refs. [25[ [3T] 
show that inclusion of a small number of couplings is al¬ 
ready sufficient for calculations in 0(V)-models. There¬ 
fore, we estimate that systematic uncertainty of our re¬ 
sults due to the truncation of 1 /( 0 ^) is comparatively 
small. 

Once we have determined the values of the critical ex¬ 
ponents, we can calculate the normalization constants 
To and Hq. The logarithms of these constants appear 
in additive terms in logM at iL = 0 and at T = Tq 
respectively, Eq. ( p!8| ). We obtain 

To = 0.0046874(8) MeV^/^ 

Ho = 13.837(13)MeV®/2 . 


Estimating the error of Tq, we perform calculations with 
the additional small explicit symmetry-breaking field: 
H = 10“^^ and 10“^^ MeV^^^. In analogy, for Hq we 
repeat our calculations in the presence of non-vanishing 
but small (T - Tc) = 1 x lO-^^ and 2 x lO-^^ MeV^/^ 
In our calculations in finite volumes an additional 
system-specific scale L should be removed in order to 
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determine finite-size scaling functions m- A possible 
normalization choice is given by 

^ = Lo\tr. (31) 

And we find 


Lo = 62.206(1) fm. (32) 

The error is estimated using the same explicit symmetry¬ 
breaking fields as in the case of error estimation for Tq. 


B. Scaling functions 

We calculate the order parameter over a wide range of 
values for T and in the presence of some small symmetry¬ 
breaking H. Our results are presented in the left part of 
the Fig. [2 If we consider a fixed temperature interval, we 
observe that for smaller values of H the order parameter 
in the vicinity of Tc decreases more rapidly and asymp¬ 
totically approaches zero already at very small tempera¬ 
ture t. It means that for small symmetry-breaking fields 
the magnetization behaves approximately as for = 0 
if we are sufficiently far away from the critical tempera¬ 
ture. However, with increasing H the crossover character 
becomes more distinct. 

We rescale the data for M using the critical exponents 
and normalization constants determined in the previous 
subsection. Our results are shown in the right part of 
Fig. [3 The rescaled data falls perfectly into one line, 
i.e., we observe ideal scaling behavior. 

For small H the scaling corrections are negligible. 
Therefore, we use the data for the smallest H to deter¬ 
mine the scaling function for the order parameter. Our 
result for /m(^) is presented in Fig. In this figure 
we also plot the order parameter scaling function for the 
0(4)-model which we obtain, using the same formalism 
as in the 0(2) case (see also [24]). We observe that these 
two functions are very similar. This similarity is a reason 
why many investigations of scaling properties of lattice 
QCD, where the lattice data is fitted to the scaling func¬ 
tions of the order parameter, led to ambiguous results 

(snaunj. 

As a check of our results, we compare our scaling func¬ 
tion in the infinite-volume limit with corresponding func¬ 
tion obtained using lattice spin simulations of 0(2)-model 
(El, Fig. i Though the calculations in m already in¬ 
clude the non-vanishing anomalous dimension we see 
an almost perfect agreement with our result. Therefore, 
we conclude that our formalism provides reasonable data. 
We infer that the RG-approach is a very appropriate tool 
for the determination of the scaling behavior. 


V. FINITE-VOLUME SCALING 

We determine the finite-volume scaling function for the 
order parameter. A similar investigation of the finite- 
size scaling in the 0(4)-model was already provided in 


m- We calculate the order parameter M as a func¬ 
tion of the symmetry-breaking field h for different finite- 
volume sizes, L =10-300 fm exactly at T = Tq. Our 
results are shown in the double-logarithmic representa¬ 
tion in the left part of Fig. In this plot we can dis¬ 
tinguish two different regions in the behavior of M: For 
very large symmetry-breaking fields we observe the same 
quantitative behavior for all volumes which we have con¬ 
sidered. The slope of the curves here is very close to 
I/(i, i.e, we observe approximately the same power law 
as in the limit L ^ oc, Eq. (18). This behavior appears 
because of the large masses of the fluctuations. In this 
situation, the correlation length ^ is much smaller than 
the extent of the system and therefore, finite L does not 
influence the critical behavior. If however, the external 
symmetry-breaking field h becomes smaller, the mass of 
fluctuations decreases and the correlation length grows. 
At some point ^ becomes so large that its value becomes 
comparable with L. Therefore, the infinite-volume scal¬ 
ing behavior becomes affected by the extent of the system 
and we observe the finite-size scaling region. 


We rescale our results and consider the rescaled order 
parameter as a function of the dimensionless scal¬ 
ing variable ^ right part of Fig. We observe that 

curves for different volumes fall almost perfectly into one 
line. However, the agreement becomes worse for decreas¬ 
ing L. This fact is explained by the presence of non- 
universal finite-size corrections. Since these corrections 
scale with the system extent as they are negligible 
for very large volumes. On account of this, we deter¬ 
mine the finite-size scaling function by using the data 
for the largest volume we have calculated (L = 300fm). 
In Fig. we present our 0(2)-finite-volume scaling func¬ 
tion together with corresponding scaling function for the 
0(4)-model. 

From the perspective of the analysis of lattice QCD re¬ 
sults, an interesting feature of the finite-size scaling func¬ 
tion is the region where the universal finite-size scaling 
behavior appears: we need to know it, in order to decide, 
where simulations need to be done in parameter space 
in order to use finite- or infinite-size scaling behavior in 
the analysis of lattice results. Also, the finite size scaling 
regions for the 0(2)- and 0(4)- models can be different 
and, probably, can be used to distinguish 0(2) and 0(4) 
universality classes. Since this regime in scaling behav¬ 
ior arises if the correlation length is in the order of the 
system extent, it is self-evident to use the universal di¬ 
mensionless combination to describe the change in 
the scaling behavior. We use the longitudinal part of the 
correlation length = I/rUcr as a measure for As also 
shown in the Fig.[^ we determine the value of this quan¬ 
tity, at the point where the finite-size scaling appears. It 
is given by 


^ = 0.395(5) for 0(2), 

/ ( 33 ) 

= 0.372(2) for 0(4). 

Lj 
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Figure 1. In this figure we present our results for the order parameter in 0(2)-model calculated in the presence of the very 
small symmetry-breaking field. In the left part we show the unrescaled, in the right part the rescaled data. The rescaled data 
fall onto one line, i.e., they exhibit the ideal universal scaling behavior. 




Figure 2. Our results for 0(2)- and 0(4)-scaling functions 
/m(^)- Both these functions behave in a very similar way, 
whereas in the vicinity of Tc the scaling function for the 0(4)- 
model decreases a bit faster with increasing than that for 
the 0(2)-model. 


Figure 3. Scaling function of the order parameter M for three 
dimensional 0(2)-model calculated within LPA using critical 
exponents from Tab. is compared with results from a lattice 
spin simulation which already include anomalous dimension 
m- We observe an excellent agreement. 


In order to determine these values, we find a linear 
parametrization for our data in a double-logarithmic rep¬ 
resentation for large and small h. We use the value of 
h at the point of intersection of these two lines as the 
point where the influence of L becomes dominant. To 
estimate the errors we determine the finite-size scaling 
region using our results for somewhat smaller volumes: 
L = 100, 200 fm. 

While the description of the finite-size scaling regions 
by the universal value of is self-evident, it is inappli¬ 
cable for practical purposes since the correlation length 
is difficult to measure in lattice simulations. Therefore 
we look for an alternative description. Since the corre¬ 
lation length of fluctuations is bounded from above by 
the inverse pion mass the finite-size scaling should 
appear if the wave-length of the pion is of the order of 


the system extent: 

— (34) 

Therefore, we can use the universal dimensionless quan¬ 
tity for the description of the point where the finite- 
size scaling occurs. Since the pion mass and the system 
extent are used as inputs in lattice simulations, this way 
to estimate the finite-size scaling region is easier to apply 
in lattice QCD. Using our approach we find: 

(m^L) = 2.01(1) for 0(2), 

= 2.12(2) for 0(4). 

Here we use the same calculation and error estimation 
technique as for f,L/L. 

The values of and rriT^L are very similar for 0(2)- 
and 0(4)-models. Therefore, using the finite-size scaling 





















Figure 4. Our results for the order parameter M in different finite volumes L = 10 — 300 fm for three-dimensional 0(2)-model 
calculated within LPA at z = 0. In the left part we present the unrescaled, in the right part the rescaled data. We see two 
different scaling regions: For large h, we see the scaling behavior valid for the limit L ^ oc, for small h we observe finite-size 
scaling. In the rescaled plot we observe that curves for different volumes fall almost perfectly onto one line. For small volumes 
we observe some deviations from the universal finite-size scaling behavior due to the non-universal scaling corrections. However, 
they are negligible for very large volumes. 





Figure 5. Finite-size scaling functions for the order parame¬ 
ter M for three dimensional 0(2)- and 0(4)-models exactly at 
the critical point t = 0 are plotted in the double-logarithmic 
representation as functions of the dimensionless scaling vari¬ 
able . The finite-size scaling regions are labeled by the 

corresponding universal values of ^rfL. For 0(2) this value 
is given by = 0.395(5), for 0(4) by = 0.372(2). 


regions to determine the nature of the chiral transition in 
lattice QCD seems to be difficult. Never the less, these 
results can be used in the scaling analysis of the lattice 
QCD data in order to clarify in which scaling region some 
particular data set should be located. Assuming 0(2) or 
0(4) scaling behavior for lattice QCD, we should fit the 
data to the finite-size scaling functions if we are in the 
region < 2 . On the other hand, for > 2 we 

can expect only infinite-size scaling. 


VI. BINDER CUMULANTS 

Next, we calculate as a function for the tempera¬ 
ture t. Since the Binder cumulant is defined in our ap¬ 
proach for finite volumes, we cannot use zero symmetry¬ 
breaking. However, we can still employ an almost vanish¬ 
ing symmetry-breaking field H = 10“^^ Mev^^^. We use 
very different volumes, L = 10 — 5000 fm. In Fig. [^we 
present our 0 ( 2 )-results for five largest volumes consid¬ 
ered. We have checked the correct asymptotic behavior 
for small T, ^4 ^ I, for all volume sizes. As expected, for 
large temperatures approaches 2 for the 0 ( 2 )-model. 
We observe the slope of B 4 decreases with decreasing L. 
For smallest volumes we have considered the correspond¬ 
ing limit can be achieved only at the temperatures far be¬ 
yond critical one. In this figure, we can see that graphs 
for different volumes cross at nearly one and the same 
point, close to the critical temperature. This crossing 
point corresponds to the universal value of the Binder 
cumulant in the thermodynamic limit. However, since 
we use large, but still finite volumes, our results include 
some finite-size corrections. We exclude them according 
to a method described in daini l33] . Our calculation 
provide following results: 

^4 = 1.2491(39) for the 0(2)-model , 

^4 = 1.0836(10) for the 0(4)-model , 

where errors are estimated using different spatial extends 
L G {1000 fm, 2000 fm, 3000 fm, 4000 fm, 5000 fm}. 

Our value for the 0(2)-models is very close to the 
value obtained in [33]: B/^ = 1.242(2). In the case 
of 0(4), we observe a somewhat larger deviation from 
the value determined using spin-model lattice simulations 
m- B/^ = 1.092(3). In both cases the discrepancy be¬ 
tween our RG-results and results from Monte Carlo sim¬ 
ulations is smaller than 1%. This deviation is explained 
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Figure 6 . Results for the Binder cumulant for the 0(2)- 
model at almost vanishing symmetry-breaking field {H — 
10"^® MeV®/2) as a function of temperature. The results 
for the largest volumes we have investigated are shown: L G 
{1000 fm, 2000 fm, 3000 fm, 4000 fm, 5000 fm}. We observe 
the correct asymptotic behavior at very large temperatures, 
B 4 ^ 2. The behavior in the limit T Tc, R 4 ^ 1, has also 
been checked. Considered on the scale of this plot, it seems 
that graphs for different volumes cross at T = Tc. 


by the fact that we neglect anomalous dimension 77. How¬ 
ever and as a matter of fact, the discrepancy is very small. 
Therefore, once again we can infer that the influence of 
finite 77 on the scaling behavior is almost negligible. 

We investigate the influence of the finite volume on 
the value of B 4 at the critical temperature. In Fig. 
we present our results for the 0(2)-model. We include 
results from data sets for all volumes we have consid¬ 
ered. We observe that for volumes which are typical for 
lattice QCD simulations 10 — 30 fm), the finite-size 
corrections to the universal value of B 4 are in the order 
of 3 — 8 % for both models. We have also checked that the 
finite-size correction is very well described by the leading 
order expansion, Eq. (27). The corresponding fit for the 
0(2)-model is also presented in Fig. (for uj we use the 
value uj = 0.6712 obtained in ED)- 

In addition, we investigate the behavior of the uni¬ 
versal value of B 4 {Tc) as a function of the symmetry¬ 
breaking field h. In Fig.j^we present our results for some 
selected volumes in rescaled form. For small h we observe 
the same finite-size effects as discussed in Fig. We do 
not observe any noticeable finite-mass corrections. For all 
volumes we find the value of B/^iTc) to remain constant 
for small h. However, at some point the value of the 
Binder cumulant decreases very fast and approaches 1. 
These two regimes are finite-size and infinite-size scaling 
regions of H 4 (Tc'): For small symmetry-breaking fields 
the correlation length is in the order of the system ex¬ 
tent L and the finite-volume effects become dominant. 
However, if we consider very large symmetry-breaking 
fields, the correlation length decreases, the finite-size ef¬ 
fects become less pronounced, and the system behaves 
approximately in the same way as in infinitely large vol¬ 
ume. With increasing h we are moving away from the 


Figure 7. Here we show the values of H 4 at — 0 for 0(2)- 
model for all volumes we have considered, L = 10 — 5000 fm. 
We observe that for volumes of some 10 fm the non-universal 
finite-size corrections to the value of the Binder cumulant at 
T = Tc become noticeable. We also provide fits of our results 
to the expansion given in Eq. (27). This form agrees very well 
with our data points. 





Figure 8. The rescaled results for B^iTc) as a function of the 
symmetry-breaking field for 0(2)-model. We observe that 
the data points do not fall onto one line. We can clearly see 
that deviations are mostly caused by non-universal finite-size 
corrections. We also present in this plot the value of ^l/L 
which correspond to the finite-size scaling region of the order 
parameter. We observe that this value describes the finite-size 
scaling region of H 4 very well. 


critical point and the system becomes more and more 
orderly. In the orderly phase, however, the Binder cumu¬ 
lant should approach the limit 1 . 

In Fig. we also plot our results for the values oi^L/L 
at the onset of the finite-size scaling regions. We find 
that our results from Fig. |^are in very good agreement 
with the finite-size scaling regions which we observe for 
the Binder cumulant. 

In Fig. we present our results for 0(2)- and 0(4)- 
models together in one plot. We observe that even in 
the presence of the finite-volume corrections, the regions 
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Figure 9. Here we present the rescaled data as in the Fig. 
but now for the 0(2)- (red points) and for the 0(4)-model 
(blue squares). For both models points with the largest val¬ 
ues of Ba correspond to calculations with L — 10 fm and 
points with the smallest values of Ba to calculations with 
L — 5000 fm. We observe that even in the presence of non- 
universal finite-size corrections we can clearly distinguish be¬ 
tween 0(2)- and 0(4)-models. 

where we can expect to measure Ba{Tc) in finite-volume 
lattice simulations do not overlap for 0(2)- and 0(4)- 
models. Of course, the smallest volume we have consid¬ 
ered is L = 10 fm, and finite-size corrections should be 
larger for even smaller volumes. However, the typical 
lattice QCD simulations are performed in volumes with 
L ^ 10 fm. Also we cannot exclude a possibility that 
non-universal finite-size corrections in the lattice QCD 
simulations are even larger than in our approach. How¬ 
ever, we see in the Fig.|^that the gap between 0(2) and 
0(4) results is relatively large. It is at least of the order of 
the finite-size corrections we expect for L = 10 fm. Thus, 
we conclude that the universal value of the Binder cumu- 
lant is a very reasonable tool to distinguish 0(2) and 0(4) 
universality classes even in the presence of finite-size cor¬ 
rections to scaling. Therefore, Ba{Tc) seems to be a very 
promising candidate for a criterion for determination of 
the universality class of the chiral transition in Nf = 2 
OT Nf = (2 + 1) lattice QCD. In order to use this re¬ 
sult, lattice QCD simulations should be performed in the 
finite-size scaling region, i.e., at the values of m-j^L < 2. 


VII. CONCLUSIONS 

We have investigated the critical behavior in the con¬ 
tinuous 0-model with 0(2)-symmetry in both infinite 
and finite volumes. In addition we have performed the 
same analysis for the 0(4) case as it was already done in 
[T^[24]. For these purposes we have applied functional 
RG approach in the local potential approximation. We 
have investigated regions where the finite-size scaling be¬ 
comes dominant. We have also studied the behavior of 


the critical fluctuations in the vicinity of the critical point 
by means of the 4^^-order Binder cumulant B^- For this 
purpose we have derived an expression for the Binder cu¬ 
mulant in the context of FRG. It allows us to investigate 
higher-order fluctuations and to better understand the 
universal and limiting behavior of B^ from a theoretical 
point of view. Though in the present work we consider 
only 0(2)- and 0(4)-models, our theoretical calculations 
can be applied to any 0(V)-symmetric model. We can 
also extend this method to Binder cumulants of higher 
order. In the case of two or (2 + l)-flavors, the lattice 
QCD simulations cannot reach exactly the critical point 
because of finite rriq and finite volumes. Therefore, ad¬ 
ditional analysis of the Binder cumulant as a function 
of the symmetry-breaking field and as a function of the 
system extent was needed in order to apply B^ for the 
analysis of lattice QCD data. 

We have considered the case of infinitely large volume 
in order to determine universal critical exponents valid in 
LPA. Our results are in a perfect agreement with results 
from [HI]. This fact implies that possible systematic un¬ 
certainty caused by the specific truncation for the order 
of the potential which we used is very small. However, 
we still cannot exclude uncertainties arising from neglect¬ 
ing the higher order kinetic terms in the ansatz for the 
scale-dependent effective action. 

We have also calculated universal scaling functions for 
the order parameter valid in infinitely large volumes. Our 
findings are in a very good agreement with results from 
where the authors have used lattice simulations for 
the three-dimensional 0(2) spin-model. The tiny devia¬ 
tions we observed arise basically because of the anoma¬ 
lous dimension r] which we have neglected in our calcu¬ 
lations. 

In our finite-volume calculations we were able to de¬ 
termine the finite-size scaling function for the order pa¬ 
rameter for 0(2)-model exactly at the critical tempera¬ 
ture. We have also described the regions where the finite- 
volume effects become dominant for 0(2)- and 0(4)- 
models using the universal values of ^l/L = 
and vrij^L. We have found that finite-size scaling re¬ 
gions for these two models are similar. Therefore, the 
difference in the finite-size scaling regions for 0(2)- and 
0(4)-models can probably not be used in order to de¬ 
termine the universality class of the chiral transition in 
Nj = 2 or Nj = (2 + 1) lattice QCD. Never the less, the 
finite-size scaling regions we have explicitly determined 
in this work are still useful for scaling analysis of lattice 
QCD data: We can use these results in order to decide 
whether a particular set of the simulation data should 
exhibit infinite-size or finite-size scaling behavior. The 
value of which separates these two regimes is given 
for both models by 2. 

We have investigated the behavior of critical fluctua¬ 
tions close to the critical point by means of the Binder 
cumulant of the 4^^ order, B^. Our numerical calcula¬ 
tions have confirmed our theoretical predictions about 
the limiting behavior of B^ and have provided the uni- 
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versal values of the Binder cumulant, exactly at the criti¬ 
cal temperature at the limit of very large volumes. These 
new FRG results are in a very good agreement with spin- 
model lattice simulations [2Ql l33] . 

Furthermore, we have investigated the influence of a 
finite symmetry breaking and finite volumes on the be¬ 
havior of B/^(Tc). We have found that in our calculations 
the finite-mass corrections are small in comparison to the 
corrections caused by the finite volumes. We have also 
shown that non-universal finite-size corrections can be 
described very well by taking only the leading order cor¬ 
rections, i.e., by contributions associated with the first ir¬ 
relevant operator in the RG-flow into account. We have 
seen that for both models such corrections are smaller 
than 8 % for volumes with L = 10 fm, which are typi¬ 
cal for lattice QCD simulations. This observation can be 
used in the analysis of Nf = 2 or Nf = (2 + 1 ) lattice 
QCD results: In Fig. we have illustrated that even in 
the presence of finite-size corrections arising in our cal¬ 
culations for 1/ = 10 fm, 0(2)- and 0(4)-models can still 
be distinguished in an unambiguous manner. Therefore, 
if we assume lattice QCD to fall into either the 0(2) or 
the 0(4) universality class, then the intervals into which 
the values for B 4 ^{Tc) measured in lattice QCD simula¬ 
tions are expected to fall are clearly different for 0 ( 2 ) 
and 0(4) universality classes. So, even for simulations 
with L ^ 10 fm we can use the universal value of the 
Binder cumulant exactly at T = Tc to determine the 
universality class of the chiral transition in Nf = 2 or 
Nf = (2 + 1) lattice QCD. However, in order to apply 
this method, lattice results should be in the finite-size 
scaling region. Also the non-universal finite-size correc¬ 
tions in lattice QCD can be potentially larger than in our 
approach. 

In conclusion, we have performed a new finite-size scal¬ 
ing analysis of the critical behavior in 0(2)- and 0(4)- 
models. Thereby we have answered the open question 
whether the finite-size scaling regions for 0(2)- and 0(4)- 
models differ or not. The difference we have observed is 
too small to be used in the scaling analysis of lattice 
QCD. In our investigation of the Binder cumulant in the 
context of the FRG, we have found that B 4 {T = Tc) 
seems to be an appropriate tool to determine the nature 
of the chiral transition in Nf = 2 or Nf = (2 + 1) lattice 
QCD simulations. We hope that these results will con¬ 
tribute fruitfully to the scaling analysis of lattice QCD. 


Appendix A: Binder cumulants from the effective 
action 


an n-point correlation function is defined as 


{pi 


113”^ 


(A 2 ) 


where Z is the generating functional and p = {cr, tt} is 
a generalized field, {p'^) contains connected and discon¬ 
nected parts. The disconnected part is given by a sum 
of products of m-point correlation functions with m < n. 
Thereby, all possible combinations with ^rrii = n ap- 

i 

pear and are multiplied with appropriate combinatorial 
factors. The connected part of a n-point correlation func¬ 
tion is given by a n^^-derivative of the generating func¬ 
tional for connected diagrams W = \ogZ: 

1 d^W 

(p")conn. = 

The first derivative of W is given by the expectation value 
of p. In following we call this quantity the classical field. 

Pel = (p) =Mp = . (A4) 


The second derivative of W is coupled to susceptibility 


d‘^W xp 1 
dHj V Vml ■ 


(A5) 


This quantity is also the dressed propagator and is 
connected to the inverse second derivative of the effec¬ 
tive action F with respect to the classical field. In our 
calculations T = VU and this statement takes the form 


d‘^W o 

_ — V^D =V 

dHp,dH^ ~ 


/ d^U \ 
\dpi^c\dp2,c\^ 


^P2 


(A 6 ) 


where Dp^p^ is a 2 x 2 matrix. 

All higher derivatives of W can be calculated using 
iterative application of the operator 


d 


E 


dp4:,c\ d 

^P4,c\ 


n 


d^U 

dp3,c\dp4:,cl 


) 


-1 


d 

dpA,c\ 


(A7) 


on the Eq. (A 6 ). 

We represent our results for the 0(2)-model in terms of 
Feynman diagrams defined as follows: A general dressed 
n-point vertex in our theory is given by 


According to Eq. (23), the 4^^-order Binder Cumulant 
for the 0 ( 2 )-model is given by 


{(j^) + (tt"^) + 2(cr^7r^) 
(cr^)^ + (Tr^)^ + 2(cr2)(7r2) 


(Al) 


Correlations appearing in this expression can be calcu¬ 
lated as follows: Eor a model with no spatial dependence 


-E 


d^U 

^01,cl • • • ^0n,cl 


d^r 

^01,cl • • • ^0n,cl 


(AS) 


and is denoted by empty circle. A general dressed static 
propagator is 


Vml 


PlP2 


(A9) 
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We denote it by line with a full circle. For a we use 
continuous and for tt dashed line. The non-vanishing 
expectation value of a is represented by 

(cr) = M = . (AlO) 

The correlation functions calculated using FRG in LPA, 
which are relevant for the of the 0(2)-symmetric 
model, are given by 


(All) 


(a^) = 





+ 4 



(A12) 


From these expressions we see that in the thermody¬ 
namic limit and in the phase with broken symmetry the 
Binder cumulant approaches 1 in the leading order. In 
the limit of restored symmetry, the order parameter M 
vanishes. Therefore, we have to look at next-to-leading 
order terms in XjV. We also have to keep in mind that 
for T ^ Tc masses of tt and a become equal. So, we 
obtain B 4 = 2 : 

In the case of 0(4)-model calculations are very simi¬ 
lar but now we have three different pion fields. However, 
they are completely equivalent. Therefore, all correlation 
functions with no mixing of different pion fields will pro¬ 
vide exactly the same results and can be calculated using 
the same expressions as for 0(2)-model. 

In correlations with mixing of pious, two different pion 
fields are involved. However any combination of different 
components will lead to one and the same result 

(ttM) = (Trlnl) , (A16) 

with i ^ j. Using this, we get the following expression 
for the Binder cumulant of the 0(4)-model: 


=. 


(A13) 


{a‘^) + 3(7r^) + 6(7r^(7^) + 6(7ri7r|) 
(a2)2+9(7r2)2 + 6((T2)(7r2) 


(A17) 


(tt") 



+ 3 ..CKD^ 



(A14) 






The only correlation in the expression above, which we 
have to evaluate in addition, is 


2 2 


) 





(AI8) 


+ 2 




Using this result we have found that limiting behavior 
of B 4 for the 0(4)-model in the thermodynamic limit is 
given by I for the low-temperature and by 3/2 for the 
high-temperature phase. 
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